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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2595 

Design  of  two-dimensional  channels  wits  prescribed 

VELOCITY  DISTRIBUTIONS  ALONG  THE  CHANNEL  WALLS 
II  -  SOLUTION  BY  GREEN'S  FUNCTION 
By  John  D.  Stanitz 


SUMMARY 

Methods . of  solution  by  Green's  function  are  developed  for  the  design 
of  two-dimensional  unbranched  channels  with  prescribed  velocities  as  a 
function  of  arc  length  along  the  channel  walls.  The  methods  apply  to 
incompressible  and  linearized  compressible,  nonviscous  irrotational  flow. 
One  numerical  example  is  presented  for  an  accelerating  elbow  with 
linearized  compressible  flow.  The  elbow  shape  obtained  from  the  solution 
by  Green's  function  is  the  same  as  that  obtained  from  a  solution  by 
relaxation  methods  for  the  same  prescribed  conditions.  The  time  required 
for  the  calculations  is  considerably  less  for  solutions  by  Green's  * 

function.  ■ 


INTRODUCTION 

,  In  this  report  a  general  method  of  design  is  developed  for  two- 
dimensional,  compressible  or  incompressible,  nonviscous  irrotational 
flow  in  unbranched  channels  with  prescribed  velocities  as  a  function  of 
arc  , length  along  the  channel  walls.  The  design  of  channels  with  pre¬ 
scribed  velocities  is  important  because:  (l)  boundary- layer  separation 
losses  can  be  avoided  by  prescribed  velocities  that  do  not  decelerate 
rapidly  enough  to  cause  separation,  (2)  shock  losses  in  compressible 
flow  and  cavitation  in  incompressible  flow  cap  be  avoided  by  prescribed 
velocities  that  do  not  exceed  certain  maximum  values  dictated  by  these 
phenomena,  and  (3)  for  compressible  flow  the  desired  flow  rate  can  be 
assured  by  prescribed  velocities  that  do  not  result  in  "choke  flow" 
conditions. 

i  _  •  . 

In  Part  I  of  this  report  (reference  1),  solutions/  were  obtained  by 
relaxation  methods.  This  method  of  solution  results  in  complete  infor¬ 
mation  concerning  the  distribution  of  flow  conditions  throughout  the 
channel  and  can  be  used  to  obtain  solutions  for  incompressible  flow  anti 
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for  two  types  of  compressible  flow:  the  general  type  with  arbitrary 
value  for  the  ratio  of  specific  heats  y  (1*4=,  for  example)  and  the 
linearized  type  with  y  equal  to  -1.0. 

In  the  present  report  solutions  are  obtained  by  Green1 s  function. 
This  method  of  solution  is  limited  to  incompressible  and  linearized 
(y  =  -1.0)  compressible  flow,  but  the  method  is  more  rapid  than  relax¬ 
ation  methods,  provided  information  within  the  channel  is  not  required. 
The  method  of  solution -is  developed  for  the  channel  walls  only  although 
the  method  can  be  extended  to  determine  the  shape  of  streamlines  within 
the  channel. 


The  design  method  reported  herein  was  developed  at  the  NACA  Lewis 
laboratory  during  1950  and  is  part  of  a  doctoral  thesis  conducted  with 
the  advice  of  Professor  Ascher  H.  Shapiro  of  the  Massachusetts  Institute 
of  Technology.  * 


METHOD  OF  SOLUTION 

The  design  method  is  developed  for  two-dimensional  channels  with 
prescribed  velocities  along  the  channel  walls.  The  prescribed  velocity 
is  arbitrary  except  that  stagnation  points  (zero  velocity)  cannot  be 
prescribed.  This  exception  limits  the  design  method  to  unbranched 
channels.  In  the  present  report  the  method  of  solution  is  by  Green  s 
function  in  conjunction  with  a  formula  derived  (elsewhere)  from  Green's 
theorem. 


Preliminary  Considerations 

Assumptions.  -  The  fluid  is  assumed  to  be  nonviscous ^and  either 
compressible  or  incompressible.  If  the  fluid  is  compressible,  the 
ratio  of  specific  heats  y  is  assumed  to  be  -1.0,  so  that  the  differ¬ 
ential  equations  describing  the  flow  are  linear.  The  flow  is  assumed  to 
be  two-dimensional  and  irrotational. 

Physical  plane.  -  The  flow  field  of  the  two-dimensional  channel  is 
considered  to  lie  in  the  physical  xy-plane  where  x  and.  y;  are 
Cartesian  coordinates  expressed  as  ratios  of  a  characteristic  length 
equal  to  the  constant  channel  width  downstream  at  infinity.  (All  symbols 
are  defined  in  appendix  A.)  < 

At  each  point  in  the  channel  the  velocity  vector  -(fig.  1)  has  a 
magnitude  Q  and  a  direction  0  where  Q  is  the  fluid  velocity 
expressed  as  the  ratio  of  a  characteristic  velocity  equal  to  the  constant 
channel  velocity  downstream  at  infinity.  For  compressible  flow,  the 
velocity  Q  is  related  to  the  velocity  q  hy 
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q  =  Qqd  (i) 

where  q  is  the  velocity  expressed  as  a  ratio  of  the  stagnation  speed 
of  sound  and  the  subscript  d  refers  to  conditions  downstream  at 
infinity. 

Stream  function  Hf .  -  If  the  condition  of  continuity  is  satisfied,, 
a  stream  function  T[r  can  be  defined  such  that  for  incompressible  flow 


where ,  from  Part  I, 


d\|r '  =s  Q  dn 


(2a) 


(2b) 


where  n  is  distance  in  the  xy-plane  measured  normal  to  the  streamline 
and  expressed  as  a  ratio  of  the  channel  width  downstream  at  infinity. 
For  linearized  compressible  flow  (y  =  -1.0) 


cl  w  -  £ 

d¥  ~  2  A\|r’ 


(2c) 


where,  from  Part  I, 


d\|r*  =  p*q*  dn  (2d) 

and  where  A\|r*  is  the  value  of  ^  *  along  the  left  channel  wall  when 
faced  in  the  direction  of  flow  if  the  value  of  r*  along  the  right  wall 
is  arbitrarily  equal  to  zero.  The  value  of  A\|r*  is  obtained  by  inte¬ 
grating  equation  (2d)  across  the  channel  at  a  position  far  downstream 
where  flow  conditions  are  uniform 


£i|r*  =  Pd*qd*  (2e) 

From  Part  I,  p*  is  related  to  the  density  p,  expressed  as  the  ratio 
of  a  characteristic  density  equal  to  the  stagnation  density,  by 
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in  which 


and 


(2h) 


(21) 


where  the  subscripts  a  and.  b  refer  to  quantities  related  to  any  two 
selected  values  of  velocity  (qa  and  q^,  respectively)  for  which  veloc¬ 
ities  the  densities  given  by  equations  (2h) -and  (2i)  are  equal  to  the 
densities  p  given  by  equations  (2f)  and  (2g) .  Also,  from  Part  I, 
q*  is  related  to  q  by 


* 


For  each  prescribed  velocity  distribution  along  the  channel  walls  there 
are  an  infinite  number  of  linearized  compressible  flow  solutions, 
depending  on  the  selected  values  of  qa  and  q-^  in  equations  (2g)  and 

(2k).  However,  for  values  of  qa  and  qb  within  the  range  of  q  pre¬ 
scribed  along  the  channel  walls  (and  therefore  everywhere  in  the  channel), 
the  solutions,  that  is,  channel  shapes,  probably  differ  only  in  small 
detail  (Part  I) . 


* 


0 
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The  values  of  qa  and  q^  might,  for  exanqple ,  he  selected  to  equal 
the  maximum  and  minimum  values  of  q  (which' values  of  q  must  occur  on 
the  channel  walls  and  are  therefore  known).  Also,  the  values  of  qa 
and  q^  might  he  selected  to  equal  the  upstream  and  downstream  veloc¬ 
ities  qu  and  q^.  In  this  case  the  upstream  and  downstream  channel 
widths  would  then  satisfy  continuity  for  a  gas  with  the  cprrect 
(arbitrary)  value  of  y  (1.4,  for  example).  If  the  upstream  and  down¬ 
stream  velocities  are  equal,  their  value  and  the  value  of  some  other 
velocity  (the  maximum  or  minimum  velocity,  for  example)  can  he  selected 
for  qa  and  q^  or,  if  desired,  qa  can  he  equal  to  q^  so  that 


and 

‘  (2m) 


Equations  (2a)  and  (2c)  define  the  stream  function  ¥  for  incom¬ 
pressible  and  linearized  compressible  flow,  respectively.  For  both 
types  of  flow  ¥  varies  from  zero  along  the  right  side  of  the  channel, 

when  faced  in  the  direction  of  flow,  to  ^  along  the  left  side  of  the 

channel. 

Velbcity  potential  $ .  -  If  the  condition  of  irrotational  fluid 
motion  is  satisfied,  a  velocity  potential  $  can  he  defined  such  that 
for  incompressible  flow 

1 

d$  =  |  dep 


(3a) 
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where,  from  Part  I, 

dSP  =  Q  ds  (3b) 

where  s  is  distance  in  the  xy-plane  measured  along  the  streamlines 
and  expressed  as  the  ratio  of  channel  width  downstream  at  infinity. 

For  linearized  compressible  flow 


- «  it  dtp’ 

~  2  Atj?5 


(3c) 


where,  from  Part  I, 

dcp*  =  q*  ds  (3d) 

Equations  (3a)  and  (3c)  define  the  velocity  potential  $  for  incom¬ 
pressible  and  linearized  compressible  flow,  respectively. 

Outline  of  design  method.  -  Solutions  for  two-dimensional  flow  are 
bo undary- value  problems.  That  is,  the  solutions  depend  on  known  con¬ 
ditions  imposed  along  the  boundaries  of  the  problem.  In  the  inverse 
problem  of  channel  design  the  geometry  of  the  channel  walls  in  the 
physical  xy-plane  is  unknown.  This  unknown  geometry  apparently  precludes 
the  possibility  of  solving  the  problem  in  the  physical  plane  and  neces¬ 
sitates  the  use  of  some  new  set  of  coordinates,  that  is,  a  transformed 
plane,  in  which  to  solve  the  problem.  These  new  coordinates  must  be 
such  that  the  geometric  boundaries  along  which  the  velocities  are  pre¬ 
scribed  are  known  in  the  transformed  plane.'  It  is  also  necessary,  for 
the  method  of  solution  employed  in  this  report,  that  the  coordinate 
system  of  the  transformed  plane  be  orthogonal  in  the  physical  plane. 

A  set  of  coordinates  that  satisfies  these  requirements  is  provided  by 
$  and  ¥,  which  are  orthogonal  in  the  physical  xy-plane  and  for  which 
the  geometric  boundaries  are  known  constant  values  of  ¥  (equal  to  0 

and  -|)  in  the  transformed  $¥-plane.  The  distribution  of  velocity  as 
a  function  of  $  along  these  boundaries  of  constant  ¥  is  known  because, 
if 

Q  =  Q(s) 
or 

'  q  =  q(s) 

is  prescribed,  equation  (3a)  or  (3c)  integrates  to  give 


0  =  $(s) 
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from  which 

Q  =  Q($) 

or 

q  =  q(<S>) 

The  technique  of  the  proposed  channel-design  method  is  therefore  to 
solve  for  the  physical  x,y- coordinates  of  the  channel  walls  in  the 
transformed  $ij/-plane  where  the  prescribed  boundary  conditions  for  the 
two-dimensional  flow  problem  are  known. 

Channel  wall  coordinates.  -  From  Part  I  the  distribution  of  channel 
wall  coordinates  x  and  y  along  the  boundaries  of  constant  ¥  equal 

to  0  and  ^  in  the  transformed  $ijr-plane  is  given  by 

(4a) 

and 


for  linearized  compressible  flow,  and  for  incompressible  flow 


(4b) 


and 


(5a) 


(5b) 


where  the  constants  of  integration  are  selected  to  give  known  (specified) 
values  of  x.  or  y  at  one  value  of  $  along  each  boundary.  Because 
q*  and  Q  are  known  functions  of  $  from  the  prescribed  velocity  as  a 
function  of  arc  length  along  the  channel  walls,  the  shape  of  the 'channel 
walls  in  the  physical  xy-plane  is  given  by  equation  (4)  or  (5)  if  6  is 
determined  as  a  function  of  $  along  the  channel  walls  equals  0  and 

tj)  .  In  this  report  the  solution  for  0  as  a  function  of  $  along  the 
channel  walls  in  the  $¥  -plane  is  obtained  by  Green's  function. 
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Solution  by  Green1 s  Function 

Continuity,  -  From  Part  I  the  continuity  equation  becomes  in  the 
transformed  <J>ur-plane 


5  loge  V  ae 

— 5$ —  +  Sf=0 

where  for  incompressible  flow 

V  =  Q 

and  for  linearized  compressible  flow 


(6a) 


(6b) 


1  + 


q* 

jl  +  q*2 


(6c) 


Irrotational  motion,  -  From  Part  I  the  equation  for  irrotational 
motion  becomes  in  the  transformed  <J>¥-plane 


a  loge  V 

5?  M  ~ 


(?) 


Integral  equation  for  0($O^1FO)  -  -  From  equations  (6a)  and  (7) 


a2e  a2e 
a$2  +  ais2 


0 


so  that  from  appendix  B  the  value  of  0  at  a  point  ($0j^0)  within, 

or  on,  the  channel  walls  in  the  transformed  ^-plane  is  given  by  the 
integral  equation 


a  log 

— w 


it 

2 


a  log  v 

— w 


d$ 
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where  the  subscripts  0  and  £  refer  to  the  channel  wall  boundaries 

2 

along  which  f  is  0  and  respectively,  and  G  is  the  Green's 
function  of  the  second  kind  for  the  channel,  which  is  an  infinite  strip 
of  width  extending  in  the  ^-direction  to  4». 


Green's  function  G.  -  The  Green's  function  of  the  second  kind  G 
for  the  infinite  channel  in  the  $¥- plane  is  given  along  the  channel 

wall  boundaries  (¥  equals  0  and  by  (appendix  C) 


G0  or  1  =  -l0Se 

O 


cosh2  ($-$0)  -  cos1 


(10) 


where  ($,W)  is  any  point  on  the  channel  wall  boundary  and 

the  point  in  the  channel  or  on  the  boundary  at  which  6  is-  to  be  deter¬ 
mined. 


Numerical  integration  for  0(^oj^o)-  -  From  equations  (9)  and  (10) 


2* 


d  loge  V 

5$ 


log. 


cosh 


(®-*0) 


sin2  ¥ 


d<3>-$  ) 


cosh 


(®-®G) 


2 

cos  ¥ 


(11) 


in  which  the  independent  variable  of  integration  has  been  changed  from 
d$  to  d($-<2?0)  so  that  the  origin,  for  purposes  of  integration,  lies 
at  $0  rather  than  $  =  0.  If  for  small  changes  in  (®.$0)  ,  that  is 

d  loge  V 

for  small  A$,  the  term  - ^ -  may  be  considered  constant  and  equal 

to  its  average  value  over  the  interval  A$,  then 


S  loge  V  A  loge  V 

— M —  =  — 


and  equation  (11)  becomes 
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2Jte($o,¥o) 


E  \ 

($-$0)=-®l 


A  loge  V 
A<2> 


($-$o)-tA$ 


loge  cosh2($-$0)  -  sin2  ¥Q  d(4>-$o)j^  - 


(A  log  V  />-®0)hA® 

L  — a#~  loge 


^  1  °e  1  "  ,cosh2($-$0)  -  cos2 


\j  d($-$o)V 


(12) 


where  the  summation  sign  is  understood  to  mean  that  the  quantity  within 
the  braces  is  summed  over  the  entire  range  of  (0-$o)  between  ±°°. 

Equation  (12)  determines  6  at  any  point  in  the  flow  field 
(channel).  For  a  point  ($o,¥q)  on  the  channel  walls  ¥Q  is  equal 

to  0  or  —  and  the  integrands  in  equation  (12)  become 
2 

2  loge  cosh  |($-$0)| 
or 

2  loge  sinh  |($-<I>0)| 

j 

so  that  equation  (12)  becomes 


($_$Q)=_oo 


'A  loge  V 
A<f> 


AI 


'A  loge  V 
A<J> 


AI 


Oj 


(13a) 


where 


M  “  J(<|)-$0)4A$  " 


(13b) 


t 
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a 

if 

2 

I*  = 

3 

if 

2 

H 

O 

11 

a 

if 

11 

O 

H 

3 

if 

■o  "  2 


J 


(13c) 


where 


a  =  ± 


log_  cosh  |($-$D)|  d|(0-$>o)| 


(13d) 


loge  sinh  |($-$)  |  d|($-$)|  (I3e) 

where  the  +  signs  apply  for  positive  values  of  ($-$0)  and  the  -  signs 
apply  for  negative  values  of  ($-$0) .  Methods  of  evaluating  a  and  p 
are  given  in  appendix  D  and  tabulated  values  are  given  for  a  wide  range 
of  |(®-®o)|  in  table  I.  Equation  (13a)  determines  0($O,YO)  &t  any 
point  on  the  channel  wall  boundaries.  Thus  from  equations  (4a)  and  (4b) 
or  (5a)  and  (5b)  the  coordinates  for  the  channel  wall  shape  in  the 
physical  xy-plane  can  be  determined. 


NUMERICAL  PROCEDURE 

The  numerical  procedure  for  the  channel  design  solution  by  Green's 
function  is  the  same,  except  for  minor  details,  for  incompressible  and 
linearized  compressible  flow.  The  stepwise  procedure  is  outlined  as 
follows: 

i 

'  (1)  For  incompressible  flow  the .velocity  Q  and  for  linearized 

compressible  flow  the  velocity  q ,  or  which  is  the  same  thing  the 
velocity  Q  and  the  constant  downstream  velocity  q^,  are  specified 
as '  functions  of  arc  length  along,  the  channel  walls  , 

Q  =  Q(s)  (14a) 
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or 

q  =  q(s)  (14b) 

where  s  is  arbitrarily  equal  to  0  at  that  point  along  one  channel 
wall  where  the  velocity  first  begins  to  vary. 

(2)  Compute  V  as  a  function  of  s  from  equations  (6b)  and  (14a) 
for  incompressible  flow  or  from  equations  (2 j) ,  (2k) ,  (6c) ,  and  (14b) 
for  linearized  compressible  flow. 

V  =  V(s)  (15) 

(3)  Compute  $  as  a  function  of  s  from  equations  (3a)  and  (3b) 

for  incompressible  flow  or  from  equations  (2e) ,  (3c) ,  (3d) ,  and  (14b) 
for  linearized  compressible  flow.  In  equation  (2e)  p^*  is  obtained 

from  equations  (2f )  to  (2i) .  For  arbitrary  distributions  of  Q  or  q 
equation  (3a)  or  (3c)  is  integrated  numerically  using,  for  example, 
Simpson's  one -third  rule.  Thus 

$  =  <J>(s)  (16) 

(4)  From  equations  (15)  and  (16)  V  and  $  are  known  functions  of 

s  so  that  * 

V  =  V($)  (1?) 

Thus  V  is  a  known  function  of  $  along  the  channel  wall  boundaries  in 

the  transformed  plane . 

(5)  If  the  prescribed  velocity  distribution  along  one  wall  is 
different  from  that  along  the  other,  the  channel  will,,  in  general,  turn 
the  flow.  This  turning  angle  A0  is  given  by  equation  (E5)  in. 
appendix  E. ‘  If  the  turning  angle  is  unsatisfactory  a  new  distribution 
of  velocity  as  a  function  of  s  (equations  (14a)  and  (14b))  is  pre¬ 
scribed  and  steps  (1)  to  (5)  repeated  until  the  desired  value  of  A0  is 
obtained.  Equation  (E5)  is  integrated  numerically  using  Simpson' s  one- 
third  rule,  for  example ,  and  equation  (17). 

(6)  The  channel  wall  boundaries  are  straight  parallel  lines  of 
constant  ¥  equal  to  0  and  £,  and  extending  to  '±=°  in  the 

u 

^-direction.  Along  these  boundaries  of  constant  ¥  a  series  of  equally 
spaced  points  are  located  at  each  of  which  the  flow  direction  Q'  and  the 
physical  x^y- coordinates  will  be  determined  by  numerical  integration. 

In  order  to  use  the  tables  of  ct  and  p  presented  in  this  report }  the 
point  spacing  must  be  an  even  multiple  of  at/24.  Thus  the  smallest 

point  spacing  ot/24  is  equal  to  l/l2  of  the  channel  width  (ot/2)  •  For 


/ 
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a  particular  prescribed  velocity  distribution  along  the  channel  walls 
the  accuracy  of  the  solution  increases,  and  so  does  the  amount  of  com¬ 
puting,  as  the  point  spacing  is  reduced.  The  error  for  a  given  point 
spacing-  depends  on  the  prescribed  velocity  distribution  and  its  order 
of  magnitude  is  given  by  the  leading  term  of  the  error  series  of  the 
formula  used  for  numerical  integration  (table  VIII,  reference  2,  for 
example) .  For  the  numerical  example  presented  in  this  report  the  point 
spacing  A<£>  was  it/12.  From  equation  (17) 

&  loge  V  (loge  -  (loge  V)t 

A$  A$>  (18) 

where  the  subscripts  $  and  $+A$  refer  to  adjacent  points  along  the 
channel  boundaries. 

(7)  The  value  of  6  at  each  point  ($0,¥Q)  on  the  channel  wall 

A  log  V 

boundaries  is  obtained  from  equation  (13a)  in  which  - ^ —  is  given 

by  equation  (18)  and  AI  is  given  by  equations  (13b),  (l3c),  and  table  I. 
Note  that  in  equation  (13a)  the  origin  has  been  moved  to  3^  by  changing 

A  loge  V 

from  $  to  ($-  <^)  .  Thus  the  value  of  - ^ —  for  a  given  value  of 

($-$Q)  varies  with  <J>0. 

(8)  The  physical  x,y- coordinates  at  each  point  on  the  channel  wall 
boundaries  are  obtained  by  the  numerical  integration  of  equations  (5a) 
and  (5b)  for  incompressible  flow,  or  equations  (4a)  and  (4b)  for  linear¬ 
ized  compressible  flow  where  A\|r*  is  given  by  equation  (2e)  .  The  con¬ 
stants  of  integration  in  equations  (4) .and  (5)  are  selected  to  give  known 
values  of  x  and  y  at  upstream  or  downstream  positions  where  flow 
conditions  can  be  considered  uniform. 


NUMERICAL  EXAMPLE 

The  channel  design  method  of  this  report  has  been  applied  to  the 
design  of  an  elbow  for  the  same  conditions  as  example  IV  of  Part  I. 

The  design  is  for  an  accelerating  elbow  with  no  local  decelerations  of 
the  prescribed  velocities  along  the  channel  walls  and  with  linearized 
compressible  flow. 

Prescribed  velocity  distribution.  -  The  prescribed  velocity  dis¬ 
tribution  along  the  channel  walls  is  given  by  q^  downstream  of  the 
elbow  and  by  Q  as  a  function  of  s  along  the  elbow  walls.  The  down¬ 
stream  velocity  qd  is  0.80176.  Along  the  inner  wall  (with  smaller 

radii)  of  the  elbow  the  arbitrarily  prescribed  velocity  Q  as  a  function 
of  arc  length  s  is  given  by 
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Q  =  0.5 


(s  <  0) 


2  3 

Q  =  i  +  -  |=  (0  <  s  <  3.0) 


Q  =  1.0 


(s  >  3.0) 


This  velocity  distribution  is  plotted  in  figure  2. 
From  equations  (l) ,  (2 j) ,  (3c),  and  (3d) 


_  *  rt  ^2^d  , 

d$  =  2  -AT*  QdS 

which  together  with  equation  (19)  integrates  to  give 


.  Jt  k2ld  (  x 

$  =  2  “Ap  (0-5  S) 


(s  <  0) 


®  =  I  "Ip  (t  +  Is  ‘  fee)  (0  -  s  -3,0) 


*  rt  k2<ld  / 

$  =  2  “5p  (-°-75+s) 


(s  >  3.0) 


where  from  equations  (2h),  (2i) ,  and  (2k)  the  constant  k2  is  equal  to 
1.36332  and  from  equations  (2e)  to  (2k)  the  constant  Ailrf  is  equal  to 
0.73782.  From  equations  (l) I,  (2j) ,  (6c),  (19),  and  (20)  the  variation 
in  loge  V  with  $  was  obtained  and  is  plotted  in  figure  3. 

The  distribution  of  velocity  as  a  function  of  arc  length  is  the 
s&me  for  both  channel  walls,  but,  as  indicated  in  figure  3,  the  dis¬ 
tribution  on  the  outer  wall  (larger  radii  in  xy-plane)  is  shifted  in  the 

positive  $ -direction  an  amount  equal  to  ^  jt  relative  to  the  distri- 

o 

bution  on  the  inner  wall.  Thus,  a  velocity  difference  exists  on  the  two 
walls  at  equal  values  of  $  in  the  interval  0  <  $  5.  3.333it,  as  shown 
in  figure  3.  The  greater  this  difference  in  velocity  or  the  greater  the 
range  in  $  over  which  velocity  differences  exist,  the  greater  is  the 
elbow  turning  angle.  For  the  prescribed  velocity  distribution  given  in 
figures  2  and  3  the  elbow  turning  angle  given  by  equation  (E5)  is 
-104.08°. 
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Results.  -  The  elbow  design  resulting  from  the  prescribed  velocities 
given  in  figures  2  and  3  is  plotted  in  figure  4.  As  indicated  in 
table  II  the  contour  of  this  elbow  is  very  nearly  the  same  as  that 
obtained  by  relaxation  methods  for  linearized  compressible  flow  with  the 
same  prescribed  conditions  (example  IV,  Part  I)  . 

The  solution  obtained  by  Green’s  function  (Part  II)  required  one 
experienced  computer  3  days  whereas  the  solution  by  relaxation  methods 
(Part  I)  required  about  10  days.  The  relaxation  solutions  provide 
additional  information,  such  as  the  distribution  of  velocity  across  the 
channel,  but  for  the  most  part  this  additional  information  is  of  second¬ 
ary  importance  and  the  design  of  channels  by  Green’s  function  is  more 
rapid  and  therefore  to  be  preferred  over  the  design  by  relaxation 
methods . 


SUMMARY  OF  RESULTS 

Methods  of  solution  by  Greek’s  function  are  developed  for  the  design 
of  two-dimensional  unbranched  channels  with  prescribed  velocities  as  a 
function  of  arc  length  along  the  channel  walls.  The  methods  apply  to 
incompressible  and  linearized-compressible,  nonviscous  irrotational  flow. 
One  numerical  example  is  presented  for  an  accelerating  elbow  with 
linearized  compressible  flow.  The  elbow  shape  obtained  from  the  solution 
by  Green’s  function  is  the  same  as  that  obtained  from  a  solution  by 
relaxation  methods  for  the  same  prescribed  conditions.  The  time  required 
for  the  calculations  was  considerably  less  for  the  solution  by  Green’s 
function. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  September  6,  1951 
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APPENDIX  A 


SYMBOLS 

The  following  symbols  are  used  in  this  report: 

,  .  ..  Bernoulli’s  numbers 

constant,  equation  (B3)  , 

Green's  function  of  the  second  kind,  equations  (B2) 
and  (10) 

integral  (a  or  p) 

coefficient,  equation  (2g) 

coefficient,  equation  (2k)‘ 

length  of  closed  boundary 

distance  in  xy-plane  measured  normal  to  direction  of  flow 
(expressed  as  ratio  of  characteristic  length  equal  to 
channel  width  downstream  at  infinity) 

velocity  (expressed  as  ratio  of  characteristic  velocity 
equal  to  constant  channel  velocity  downstream  at 
infinity) 

velocity  (expressed  as  ratio  of  stagnation  speed  of 
sound) 

velocity  used  in  linearized  compressible  flow  and  related 
to  q  by  equation  (2j) 

distance  from  any  point  in  $¥-plane  to  point  ($0,¥Q) 
at  which  logarithmic  singularity  exists 

distance  in  xy-plane  measured  along  direction  of  flow 
(expressed  as  ratio  of  characteristic  length  equal  to 
channel  width  downstream  at  infinity) 

velocity  parameter  defined  by  equations  (6b)  and  (6c)  for 
incompressible  and  linearized  compressible  flow, 
respectively 


NACA  TS  2595 


17 


V,v1}v2 

complex  functions  defined  by  equations  (C3),  (cla),  and 
(C2a) ,  respectively 

x,y 

Cartesian  coordinates  in  physical  plane  (expressed  as 
ratios  of  characteristic  length  equal  to  channel  width 
downstream  at  infinity) 

z 

complex  coordinate,  equation  (Clb) 

z 

conjugate  of  z 

a 

integral,  equation  (13d) 

P 

integral,  equation  (l3e) 

r 

ratio  of  specific  heats 

A 

finite  increment 

e 

flow  direction  in  physical  xy-plane  (measured  in  counter¬ 
clockwise  direction  from  positive  x-axis) 

a e 

channel  turning  angle,  equation  (El) 

p 

density  (expressed  as  ratio  of  stagnation  density) 

p* 

density  used  in  linearized  compressible  flow  and  related 
to  p  by  equation  (2f) 

velocity  potential  used  as  Cartesian,  coordinate  in  trans¬ 
formed  $¥-plane  and  related  to  cp  or  cp*  by  equa¬ 
tion  (3a)  or  (3c) ,  respectively 

cp  and  cp* 

velocity  potential  for  incompressible  and  linearized 
compressible  flow,  respectively,  equations  (3b)  and  (3d) 

¥ 

stream  function  used  as  Cartesian  coordinate  in  trans¬ 
formed  $¥-plane  and  related  to  \Jr  Or  ¥*  by  equa¬ 
tion  (2a)  or  (2c) ,  respectively 

¥  and  * 

stream  function  for  incompressible  and  linearized  com¬ 
pressible  flow,  respectively,  equations  (2b)  and  (2d) 

Ai|r* 

boundary  value  of  i|r*,  for  linearized  compressible  flow, 
along  left  channel  wall  when  faced  in  the  direction  of 
flow,  equation  (2e) 

18 
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Subscripts:. 

a,b 

/ 

d 

o 

XL 

{*-%) 

($_$Q)  +  A$ 
0 
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any  harmonic  function  in  $  If- plane 

quantities  related  to  any  two  selected  values  of  velocity 
(qa  and  qb,  respectively)  for  which  densities  given 
hy  equations  (2h)  and  (2i)  are  equal  to  densities  given 
hy  equations  (2f)  and  (2g) 

conditions  downstream  at  infinity 

point  in  $¥-plane  at  which  Q  is  determined 

conditions  upstream  at  infinity 

J>oint  at  ($-$Q)  on  either  channel  wall  boundary 
point  at  £($-  $o)  +  At]  on  either  channel  vail  boundary 
boundary  along  which  If  equals  0 
boundary  along  which  f  equals  ^ 


it 

2 
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APPENDIX  B 


INTEGRAL  EQUATION  FOR  0($o,¥o) 

If  the  distribution  of  the  angle  0($,V)  in  the  transformed 
W-plane  is  harmonic ,  that  is ,  satisfies  equation  (8)  within  and  on  the 

channel  walls  (¥  equals  0  and  —) ,  then  from  Green's  theorem  and 

the  theorem  of  mean  value  it  can  be  shown  that  the  value  of  6  at  a 
point  within  (or  on)  the  channel  walls  is  given  by  (refer¬ 

ence  3 ,  p.  204,  for  example) 


(Bl) 


where  the  two  integrals  on  the  right  side  of  equation  (Bl)  represent 
the  line  integral  around  the  channel  walls  in  the  counterclockwise 

direction  with  the  signs  adjusted  so  that  represents  the  inner 

o¥ 

normal  to  the  path  of  integration. 


The  function  G($>,Uf)  in  equation  (Bl)  is  of  the  form  (refer¬ 
ence  3,  p.  204). 


G(M  =  loge  i  +  o>(S>,¥) 


(B2) 


where  r  is  the  distance  frpm  any  point  ($,1?)  to  the  point  ($0>¥0) 

and  where  o)($,¥)  is  an  arbitrary  function  that  is  harmonic  within  and 
on  the  channel  walls.  (Thus  from  equation  (B2) ,  G($,¥)  is  harmonic 

within  and  on  the  channel  walls  except  at  the  point  ($0^0)  where  a 
logarithmic  singularity  exists.)  Because  the  harmonic  function 
is  arbitrary,  the  function  G($,¥)  can  be  selected  so  that  along  the 

channel  wall  boundaries  (¥  equals  0  and  — )  —  is  a  constant  c 

*  d¥ 

given  by  the  following  equation  (obtained  from  notes  presented  by 
Tamarkin  and  Feller  in  the  1941  Summer  Session  for  Advanced  Instruc¬ 
tion  and  Research  in  Mechanics  at  Brown  Univ.): 


where  l  is  the  length  of  the  path  along  which  the  line  integral  is 
taken.  For  the  path  under  consideration  l  is  infinite  and  therefore 
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G($,V)  can  be  selected  so  that  ^2  is  zero  along  the  channel  walls. 

A  function  with  this  property  is  called  a  Green*  s  function  of  the  sec 
ond  kind.  Equation  (Bl)  he  comes 


or,  combined  with  equation  (6a) 


d  log  V 

Along  the  channel  walls  - — s —  is  known  from  the  prescribed  veloc- 

S  $ 

ity  distribution  so  that,  after  the  proper  Green' s  function  G  has 
been  determined  (appendix  C),  equation  (9)  determines  the  value  of  9 
at  any  point  ($0,¥0) •  The  value  of  0($O,WO)  given  by  equation  (9) 

can  be  adjusted  by  an  arbitrary  constant  of  integration  to  give  a  spec¬ 
ified  value  of  9  at  one  point  in  the  flow  field. 
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APPENDIX  C 


GREEN’S  FUNCTION  OF  SECOND  KIND 


From  appendix  B  Green’s  function  of  the  second  kind  G  satisfies 
the  condition 

-  0 

along  the  channel  vails ,  which  are  straight  and  parallel  boundaries 
(»  equals  0  and  ~)  extending  to  ±co  in  the  $  -direction,  and 
satisfies  the  equation 


everywhere  in  the  channel  except  at  the  point  .($0'^o)  w^ere  0  has 

a  logarithmic  pole.  For  these  conditions  the  Green’s  function  G  can 
be  obtained  by  analogy  from  the  velocity  potential  for  incompressible 
flow  into  a  point  sink  at  ($0>¥0)  between  straight  parallel  bound¬ 
aries  at  TJT  equal  to  0  and  The  logarithmic  pole  for  G  at 

^  dG 

($0,ijr0)  corresponds  to  the  point  sink  and  the  condition  —  =*  0  at 

the  boundaries  corresponds  to  zero  velocity,  that  is,  no  flow  normal 
to  the  boundaries. 


,  The  velocity  potential  for  fluid  flow  with  the  boundary  conditions 
just  described  is  obtained  from  two  infinite  series  of  point  sinks  with 
the  sinks  of  each  series  spaced  it  distance  apart  in  the  ^-direction 
and  the  two  series  arranged  by  the  method  of  images  in  such  a  manner 

that  no  flow  crosses  the  boundaries,  that  is  ^  =  0.  This  arrangement 

oW 

of  point  sinks  is  shown  in  figure  5. 


The  complex  function  wj_  for  the  first  infinite  series  of  point 
sinks  is  given  by  (reference  4,  p.  112,  for  example) 

v1  =*  -loge  sinh  (z-zQ)  (Cla) 


where 


(Clb) 


Z  a  ■$  +  i*JT 
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The  complex  function  w2  for  the  second  infinite  series  of  point  sinks 
(mirror  image  of  the  first  series  in  order  to  prevent  flow  across  the 
"boundaries  ¥  equals  0  and  is  given  by  ' 

w2  =  -loge  sinh  (z-zQ)  (C2a) 

where 

z  =  $  -  i¥  (C2b) 

The  complex  function  w  for  the  combined  flow  becomes  from  equa¬ 
tions  (Cla)  to  (C2b) 

w  =  w-l  +  w2  =  -loge  sinh  [($-  $Q)  +  i  (¥-¥0j]  - 

loge  sinh  +  1  (¥+Wo)3  (C3) 

The  Green* s  function  of  the  second  kind  G  corresponds  to  the 
velocity  potential  for  the  incompressible  flow  and  is  therefore  given  by 
the  real  part  of  equation  (C3) 

G  =  -  ~  loge  jcosh2  ($-$0)  -  cos2  (¥-¥Q)J  [cosh2  ($-  $Q)  -  cos2  (¥+¥0)| 

(04) 

It 

But  along  the  channel  walls  ¥  is  equal  to  0  or  so  that 

cos2'  (M0)  =  cos2  (¥-¥0) 
and  equation  (C4)  becomes 

G0  or  1  =  _loSe  [cosh2  ($-<J>0)  -  cos2  (¥-¥0)]  (10) 

2  ‘  ' 

Equation  (10)  gives  the  Green’s  function  of  the  second  kind  along  the 
channel  walls  (straight  parallel  lines  of  constant  ¥  equal  to  0  and 

^  and  extending  to  in  the  $- direction) . 

2  .  • 


NACA  TN  2595 


23 


APPENDIX  D 


EVALUATION  OF  a  AND  0 

Several  techniques,  depending  on  the  magnitude  of  the  upper  limit 
|(®-M  ,  were  used  to  evaluate  the  integrals  a  and  0  given  by 
equations  (13d)  and  (l3e) .  Each  integral  is  treated  separately  in  this 
appendix  and  the  values  of  ($-.$0)  for  the  'upper  limit  |($-$0)j  are 
considered  positive.  For  negative  values  of  (3>-$0)  the  magnitudes 

of  I  (that  is,  of  a  or  0)  are  equal  for  corresponding  values  of 
|($-$0)|  but  opposite  in  sign.  As  a  result  the  values  of  AI  have  the 
same  sign. 


Integral  a 

Small  and  medium  values  of  ($-$0)  •  -  For  small  and  medium  values 
of  the  upper  limit  of  integration  ($-$Q)  in  equation  (13d),  that  is, 
for  0  <  (3>-$0)  <60  rt/24,  the  integral  a  is  evaluated  by  Simpson1  s 
one-third  rule  using  increments  of  ($-$  )  equal  to  jt/48. 

Large  values  of  -  For  large  values  of  ($-$0)  that  is 

for  (<£-  >  60  it/ 24,  the  integrand  in  equation  (13d)  becomes 


loge  cosh  «  ($-%)  -  loge  2 

so  that  equation  (13d)  becomes 
60*/24 


(Dl) 


a  z 


log  cosh  ($-<!>)  d($-$ »)  + 


($-%) 

60jt/24 


[(®-®Q)  -  l°Se  2]  d($-$o) 


25.809782 


4 


($-$j 


-  0.693147  -  25.398552 


■] 


*  0.411230  -  0.693147  ($-  $  )  +  i  ($>-$o)2 


(D2) 
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Equation  (D2)  gives  values  of  a  for  values  of 
greater  than  60it/24 .  Values  of  the  integral  a 
for  a  range  of  j($-$Q)|  between  0  and  100rt/24 

For  negative  values  of  ($-$0)  the  sign  of  a 

Integral  (5 

For  ($-$0)  equal  to  zero  the  integrand 

infinity  so  that  Simpson’s  one-third 
3  in  this  region  of  ($-<J>0) ,  as  was 
(l3e)  integrates  by  parts  to  give 


sinh  (®-®0)  a(®-®0) 
=  ($-4>0)  loge  sinh  ($-$Q)  - 


Small  values  of  ($-$0)  .  - 

of  equation  (l3e)  becomes  minus 
rule  cannot  "be  used  to  evaluate 
done  for  a.  However,  equation 

ioge 


equal  to  or 

are  tabulated  in  table  I 
in  increments  of  jf/24. 

is  negative. 


($-$0)  ctnh  ($>-$0)  d($-$Q) 


(D3) 


where  the  integrand  ($-$0)  ctnh  ($-$0)  on  the  right  side  of  equa¬ 
tion  (D3)  can  be  expanded  in  the  following  series  form: 


($-$o)  ctnh  (®-*0)  =  1  + 


22Bx  ($-$q)2  2% 


($-$o) 


2! 


4! 


2^  ($_$o) 

61 


28B. 


(0-®o) 


8 


210B, 


($_$o) 


10 


212B 


11 


($.$o) 


12 


8! 


10 


12 


(D4) 


where  Bj_, B3,  and  so  forth,  are  called  Bernoulli’s  numbers  (refer¬ 
ence  5,  p.  90,  for  example).  From  equations  (D3)  and  (D4) 
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P  =  (<J>-$Q)  loge  sinh  ($-^0) 


($-$D)3  ($>-%) 5  2($-$c)7 

9  +  225  "  6615  + 


($_$o)9  2  11  1382  ($-$o)13 
42,525  1,029,105  +  8,300,667,375 


(D5) 


Equation  (D5)  was  used  to  obtain  p  as  a  function  of  ($-$q)  for 
0  <  (®_®o)  <  8it/24. 


tfedium  values  of  ($-4>0)  .  -  For  medium  values  of  the  upper  limit 
of  integration  ($-  $Q)  in  equation  (l3e),  that  is,  for  8n/24  < 

<  60rt/24,  the  integral  p  is  evaluated  by  Simpson’s  one-third  rule  as 
was  done  for  a. 


Large  values  of  ($-$^) .  -  For  large  values  of  (3>-$q),  that  is, 
for  ($-$0)  >  60«/24,  the  integrand  in  equation  (l3e)  becomes 


loge  sinh  (^-$0)  *  ($-$0)  -  loge  2 

so  that  equation  (I3e)  becomes 
60jc/24 


p  s 


loge  sinh  ($-$Q)  d^-^)  + 


(D6) 


60rt/24 


[($-%)  -  logg  2]  d(®-®0) 


24.576082 


4- 


($-$>) 


0.693147  ($-$0)  -  25.398552 


3 


*»-  0.822470  -  0.693147  ($-$0)  +|  (<D-<I>0)2  (D7) 

Equation  (D7)  gives  values  of  ,p  for  values  of  (3>-$o)  equal  to  or 

greater  than  60jt/24.  Values  of  the  integral  p  are  tabulated  in  table  I 
for  a  range  of  K*-®0)i  between  0  and  100rt/24  in  increments  of  jt/24. 

For  negative  values  of  ($-$Q)  the  sign  of  P  changes. 
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APPENDIX  E 

CHANNEL  TURNING  ANGLE 

If  the  prescribed  velocity  distribution  along  one  channel  wall 
differs  from  the  distribution  along  the  other  wall,  then  in  general  the 
channel- deflects  an  amount  AS,  which'  is  the  difference  in  flow  direc¬ 
tion  far  'downstream  and  far  upstream  of  the  region  in  which  the  pre¬ 
scribed  velocity  distribution  varies.  Thus, 

A0  =  0d  -  9U  (El) 

For  large  values  of  |($-$0)|  such  as  occur  far  upstream  and  far 

downstream  of  the  region  in  which  the  prescribed  velocity  varies  along 
the  channel  walls 

cosh^  ($-$0)  »  cos^  (TJ*-¥0) 
so  that  from  equation  (10) 

G0  =  G£  =  2  [i<$-*oy|  -  loge  2  (E2) 

2 

Far  upstream  so  that 

)(®-®0)l  =  (®-®0> 

and  because  V  is  harmonic 


so  that  equation  (E2)  substituted  in  equation  (9)  gives 
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Likewise,  far  downstream  $0>  0  so  that 

|($-$0)|  =  -  (*-®0) 

and  equation  (E2)  substituted  in  equation  (9)  gives 


From  equations  (El),  (E3),  and  (E4) 


Equation  (E5)  determines  the  channel  turning  angle  AO. 


REFERENCES 

1.  Stanitz,  John  D.:  Design  of  Two-Dimensional  Channels  with  Prescribed 

Velocity  Distributions  Along  the  Channel  Walls.  I  -  Relaxation 
Solutions.  NACA  TN  2593,  1952. 

2.  Southwell,  R.  V.:  Relaxation  Methods  in  Theoretical  Physics. 

Clarendon  Press  (Oxford),  1946. 

3.  Osgood,  William  Fogg:  Functions  of  a  Complex  Variable.  G.  E. 

Stechert  &  Co.  (New  York),  1942. 

4.  Streeter,  Victor  L.:  Fluid  Dynamics.  McGraw-Hill  Book  Co.,  Inc. 

(New  York),  1948.  t 

5.  Peirce,  B.  0.:  A  Short  Table  of  Integrals.  Ginn  and  Company 

(Boston),  1929. 


28 


HACA-TN  2595 


TABLE  I  -  TABULATED  VALUES  OF  THE  INTEGRALS  a.  AMD  B  FOB  A  RANGE  OF  |(  <S-3b  )l 
[Computational  methods  given  in  appendix  D.] 


1 <*-*„>  I- 

_  ( x) 

CL 

0 

0 

1U/24) 

0.000373 

2(rt/24) 

.002970 

3(it/24) 

.009942 

4(*/24) 

.023301 

5  (jt/24) 

.044875 

6(jt/24) 

.076277 

7(*/24) 

.118897 

8(rt/24) 

. 173909 

9(  jt/24) 

.242283 

10  (jt/24) 

.324812 

ll(«/24) 

.422136 

12(  jt/24) 

.534763 

13  (jt/24) 

.663098 

14(ji/24) 

.807460 

15(n/24) 

.968096 

16  (*/24) 

1.145201 

17  (jt/24) 

1.338926 

18  (jt/24) 

1.549389 

19  (jt/24) 

1.776679 

20(jt/24) 

2.020867 

2l(  Jt/24) 

2.282008 

22(jt/24) 

2.560143 

23(jt/24) 

2.855304 

24U/24) 

3.167516 

25(rt/24) 

3.496799 

26(  Jt/24) 

3.843168 

27  (jt/24) 

4.206633 

28(  jt/24) 

4.587203 

29  (jt/24) 

4.984885 

30  (jt/24)  , 

5.399685 

31(rt/24) 

5.831607 

.  32  (jt/24) 

6.280652 

33(  Jt/24) 

6.746825 

34  (jt/24) 

7.230126 

35(ji/24) 

7.730557 

36(  Jt/24) 

8.248119 

37(ir/24) 

8.782813 

38( n/24) 

9.334640 

39  (jt/24) 

9.903600 

40(  Jt/24) 

10.489693 

41(  Jt/24) 

11.092920 

42  (jt/24) 

11.713281 

43(  Jt/24) 

12.350777 

44(  j(/24) 

13.005406 

45  (jt/24) 

13.677170 

46(jt/24) 

14.366068 

47(  Jt/24) 

15.072101 

48  ( jt/24) 

15.795268 

49  (jt/24) 

16.535571 

... 

AI  =  Aa 


A-jCt, 

(A&=*/24 


A2a 
)|(A4>=2j(/24) 


0.000373 


.002597 


.006972 


.013359 


.021574 


.031402 


.042620 


.055012 


.068374 


.082529 


.097324 


.128335 


.144362 


.160636 


193725 


.210463 


.244188 


.278135 


.295161 


.312212 


.329283 


.246369 


.363465 


.380570 


.397682 


.414800 


.431922 


.  449045 


.466173 


.483301 


.500431 


.517562 


.534694 


.551827 


.568960 


.586093 


.620361 


.637496 


.671764 


.688898 


.723167 


.740303 


.757436 


0.002970 


.020331 


.052976 


.097632 


.150903 


.209951 


.272697 


.337741 


.404188 


.471478 


.539276 


,607373 


.675652 


.744035 


.812482 


.880967 


.949474 


1.017993 


1.086521 


1.155053 


1.223588 


1.292125 


1.429200 


1.497739 


p(D 

0 

-0.396937 

-.611660 

-.756406 

-.854451 

-.916486 

-.948530 

-.954346 

-.936438 

- . 896542 

- . 835896 

-.755399 

-.655715 

-.537338 

-.400640 

-.245901 

-.073335 

. 116897 

.324671 

.549892 

.792488 

1.052403 

1.329594 

1.624029 

1.935683 

2.264536 

2.610573 

2.973783 

3.354158 

3.751689 

4.166373 

4.598205 

5.047181 

5.513301 

5.996561 

6.496961 

7.014499 

7.549175 

8.100987 

8.669936 

9.256021 

9.859241 

10.479597 

11.117089 

11.771715 

12.443477 

13.132373 

13.838405 

14.561571 

15.301873 

AI  »  Ap 


AiP 

(A$>=jt/24) 

A20 

(A$=2  jt/24) 

-0.396937 

-0.611660 

-.214723 

-.144746 

-.242791 

-.098045 

-.062035 

-.094079 

-.032044 

-.005816 

.012092 

.017908 

.039896 

.100542 

.060646 

.080497 

-.180181 

.099684 

'.118377 

.255075 

.136698 

.154739 

.327305 

.172566 

. 190232 

.398006 

.207774 

.225221 

.467817 

.242596 

.259915 

.537106 

.277191 

.294435 

.606089 

.311654 

.328853 

.674890 

.346037 

.363210 

.743585 

.380375 

.397531 

.812215 

.414684 

'431832 

.880808 

.448976 

.466120 

.949380 

.483260  ' 

.500400 

1.017938 

.517538 

.534676 

1.086488 

.551812 

.568949 

1.155034 

.586085 

.603220 

1.223576 

.620356 

.637492 

1.292118 

.654626 

.671762 

1.360658 

.688896 

. 706032 

1  1.429198 

.723166 

.740302 

1.497738 

.757436 

^For  negative,  values  of  the  signs  of  a  and  0  change,  "but  the  signs 

of  Aa  and  A0  remain  unchanged. 
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TABLE  I  -  TABULATED  VALUES  QF  THE  INTEGRALS  a  AND  p  FOR  A  RANGE  OF  |($-$0)|  -  Concluded 
[Computational  methods  given  in  appendix 


>> 

50( fl/24) 

17.293007 

5l(n/24) 

18.067579 

52(it/24) 

18.859285 

53(«/24) 

19.668125 

54( n/24) 

20.494101 

55(  it/24) 

21.337211 

56(it/24) 

22.197456 

57(it/24) 

23.074835 

58( n/24) 

23.969349 

59(it/24) 

24.880998 

60(  it/24) 

25.809782 

6l(  it/24) 

26:755694 

62(  it/24) 

27.718750 

63(it/24) 

28.698940 

64(  it/24) 

29.696257 

65(  it/24) 

30.710717 

66( rt/24) 

31.742311 

67  ( it/24) 

32.791033 

68(it/24) 

33.856896 

69  ( it/24) 

34.939895 

70{  it/24) 

36.040029 

71(it/24) 

37.157289 

72(it/24) 

38.291692 

73(  it/24) 

39.443230 

74(it/24) 

40.611893 

75(  it/24) 

41.797701 

76(it/24) 

43.000643 

77 (rt/24) 

44.220710 

78(it/24) 

45.457922 

79  (it/24) 

46.712269 

80(  it/24) 

47.983750 

Ql(  it/24) 

49.272356 

82(  it/24) 

50.578107 

83(it/24) 

51.900992 

84 (rt/24) 

53.241002 

85(«/24) 

54.598158 

86 (rt/24) 

55.972447 

87(it/24) 

57.363861 

88  ( it/24) 

58.772421 

89  (it/24) 

60.198115 

90( rt/24) 

61.640933 

91(  it/24) 

63.100897 

92  (it/24) 

64.577995 

93(  it/24) 

66.072228 

94(  it/24) 

67.583585 

95(it/24) 

69.112088 

96(rt/24) 

70.657725 

97(it/24) 

72.220486 

98(  it/24) 

73.800393 

99(it/24) 

75.397435 

100(rt/24)(2) 

77.011599 

AI  »  Ax 


A^a 
(A($=rt/24)| 


.774572 


.791706 


.808840 


.825976 


.843110 


.860245 


.877379 


.894514 


.911649 


.928784 


.945912 


,963t)56 


.980190 


.997317 


1.014460 


1.031594 


1.048722 


1.065863 


1.082999 


1.100134 


1.117260 


1.134403 


1.151538 


1.168663 


1.185808 


1.202942 


1.220067 


1.237212 


1.254347 


1.288606 


1.305751 


1.322885 


1.340010 


1.374289 


1.391414 


1.408560 


1.425694 


1.442818 


1.459964 


1.477098 


1.494233 


1.511357 


1.528503 


1.545637 


1.562761 


1.579907 


1.597042 


1.614164 


Aga 

(A4t2rt/24) 


1.566278 


1.634816 


1.703355 


1.840433 


1.908968 


1.977507 


2.046044 


2.114585 


2.183123 


2.251663 


2.320201 


2.388740 


2.457279 


2.525818 


2.662895 


2.731435 


2.799974 


2.868512 


2.937052 


3.005590 


3.074130 


3.142668 


3.211206 


p(D 

16.059309 

16.833880 

17.625585 

18.434426 

19.260401 

20.103511 

20.963755 

21.841135 

22.735649 

23.647298 

24.576082 

25.521994 

26.485050 

27.465240 

28.462557 

29.477017 

30.508611 

31.557333 

32.623196 

33.706195 

34.806329 

35.923589 

37.057992 

38.209530 

39.378193 

40.564001 

41.766943 

42.987010 

44.224222 

45.478569 

46.750050 

48.038656 

49.344407 . 

50.667292 

52.007302 

53.364458 

54.738747 

56.130161 

57.538720 

58.964415 

60.407233 

61.867197 

63.344295 

64.838528 

66.349885 

67.878388 

69.424025 

70.986786 

72.566693 

74.163735 

75.777899 

AI  = 

A0 

Ax0 

(A<k=rt/24) 

V 

(Afc=2rt/24) 

.774571 

1.566276 

.791705 

.808841 

1.634816 

.825975 

.843110 

1.703354 

.860244 

.877380 

1.771894 

.894514 

.911649 

1.840433 

.928784 

.945912 

1.908968 

.963056 

.980190 

1.977507 

.997317 

1.014460 

2.046044 

1.031594 

1.048722 

2.114585 

1.065863 

1.082999 

2.183123 

1.100134 

1.117260 

2  251663 

1.134403 

1.151538 

2.320201 

1.168663 

1.185808 

2.388740 

1.202942 

1.220067 

2.457279 

1.237212 

1.254347 

2.525818 

1.271481 

1.288606 

2.594357 

1.305751 

1.322885 

2.662895 

1.340010 

1.357156' 

2.731435 

1.374289 

1.391414 

2.799973 

1.408560 

1.425694 

2.668513 

1.442818 

1.459964 

2.937052 

1.477098 

1.494233 

3.005590 

1.511357 

1.528503 

3.074130 

1.545637 

1.562761 

3.142668 

1.579907 

1.597042 

3.211206 

1.614164 

(1) 


(2) 


For  negative  values  of  (4>-$0)  "the  signs  of  a  and 
of  Aa  and  AS  remain  unchanged. 

For  values  of  |{$-$0)|  >100(rt/24)  use  equation  (D2) 
for  B. 


0  change,  hut  the  signs 
for  a  and  equation  (D7) 
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TABLE  II  -  COMPARISON  OF  ELBOW  DESIGNS  OBTAINED  FROM  SOLUTIONS  BY  RELAXATION  METHODS  AND  BY  GREEN'S  FUNCTION 


[Linearized  compressible  flow;  prescribed  velocity  distribution  given  in  figs.  2  and  3j] 


* 

\ 

V  =  0 
(inner  wall) 

V  -  «/2 
(Outer  wall) 

Q 

<1 

Solution  by 
relaxation  methods 
(Part  I) . 

Solution  by 
Green’s  function 
(Part  II) 

Solution  by 
relaxation  methods 
(Part  I) 

Solution  by 

Green's  function 
(Part  II) 

X 

y 

e 

(deg)' 

X 

y 

e 

(deg) 

Q 

X 

y 

e 

(deg) 

X 

y 

e 

(deg) 

-22(rt/24) 

0.5000 

0.4009 

-2.466 

-0.769 

0 

-2.466 

-0.769 

0 

0.5000 

0.4009 

-2.466 

.770 

0 

-2.466 

.770 

0 

-20(rt/24) 

.5000 

.4009 

-2.241 

-.769 

.01 

-2.241 

-.769 

.01 

.5000 

.4009 

-2.241 

.770 

-.01 

-2.241 

.770 

-.01 

-18U/24) 

.5000 

.4009 

-2.016 

-.769 

.01 

-2.016 

-.769 

.01 

.5000 

.4009 

-2.016 

.770 

-.01 

-2.016 

.770 

-.01 

-16U/24) 

.5000 

.4009 

-1.791 

-.769 

.02 

-1.791 

-.769 

.02 

.5000 

.4009 

-1.791 

.770 

-.02 

-1.791 

.770 

-.02  ' 

-14(it/24) 

.5000 

.4009 

-1.566 

-.768 

.  .03 

-1.566 

-.768 

.03 

.5000 

.4009 

-1.566 

.770 

-.03 

-i.566 

.770 

-.03 

-12  ( it/24) 

.5000 

.4009 

-1.341 

-.768 

.05 

-1.341 

-.768 

.05 

.5000 

.4009 

-1.341 

.769 

-.05 

-1.341 

.770 

-.05 

-10(ir/24) 

.5000 

.4009 

-1.116 

-.768 

.08 

-1.116 

-.768 

.08 

.5000 

.4009 

-1.116 

.769 

-.08 

-1.116 

.769 

-.08 

-8(«/24) 

.5000 

.4009 

-.891 

■  -.768 

.14 

-.891 

-.768 

.14 

.5000 

.4009 

-.891 

.769 

-.13 

-.891 

.769 

-.13 

-6  ( it/24) 

.5000 

.4009 

-.666 

-.767 

.24 

-.666 

-.767 

.24 

.5000 

.4009 

-.666 

.768 

-.21 

-.666 

.768 

-.22 

-4(  it/24) 

.5000 

.4009 

-.441 

-.766 

.40 

-.441 

-.766 

.41 

.5000 

.4009 

-.441 

.767 

-.35 

-.441 

.767 

-.36 

-2  ( it/24) 

.5000 

.4009 

-.216 

-.763 

.70 

-.216 

-.763 

.74 

.5000 

.4009 

-.216 

.765 

-.56 

-.216 

.765 

-.59 

0 

.5000 

.4009 

.008 

-.760 

1.31 

.010 

-.759 

1.52 

.5000 

.4009 

.009 

.763 

-.92 

.009 

.762 

-.94 

2(it/24) 

.5079 

.4072 

.233 

•  -.752 

2.82 

.233 

-.751 

2.76 

.5000 

.4009 

.234 

.758 

-1.45 

.234 

.758 

-1.48 

-  4(  it/24) 

.5293 

.4243 

.450 

-.739 

3.88 

.449 

-.738 

3.76 

.5000 

.4009 

.459 

.751 

-2.22 

.459 

.750 

-2.24 

6(n/24) 

.5599 

.4489 

.656 

-.724 

4.00 

.656 

-.724 

3.89 

.5000 

.4009 

.684 

.740 

-3.28 

.684 

.740 

-3.29 

8(it/24) 

.5962 

.4780 

.851 

-.712 

3.17 

.850 

-.713 

13.07 

.5000 

.’4009 

.908 

.725 

-4.65 

.908 

.724 

-4.68 

10{it/24) 

.6354 

.5094 

1.033 

-.704 

1.44 

1.033 

-.705 

1.38 

.5000 

.4009 

1.132 

.703 

-6.41 

1.132 

.702 

-6.43 

12(it/24) 

.6754 

.5415 

1.205 

-.703 

-.98 

1.205 

-.704 

-1.04 

.5000 

.4009 

1.355 

.674 

-8.52 

1.355 

.673 

-8.56 

14(it/24) 

.7149 

.5732 

1.366 

-.710 

-4.00 

1.367 

-.711 

-4.04 

.5000 

.4009 

1.577 

.636 

-11.04 

1.577 

.635 

-11.06 

16(it/24) 

.7531 

.6038 

1.519 

-.725 

-7.48 

1.519 

-.727 

-7.51 

.5000 

.4009 

1.797 

.587 

-13.91 

1.797 

.586 

-13.94 

18(  it/24) 

.7894 

.6329 

1.663' 

-.749 

-11.35 

1.663 

-.750 

-11.37 

.5000 

.4009 

2.013 

.527 

-17.13 

2.013 

.526 

-17.16 

20(it/24) 

.8235 

.6602 

1.798 

-.781 

-15.52  . 

1.799 

-.783 

-15.54 

.5000 

.4009 

2.226 

.455 

-20.67 

2.226 

.453 

-20.69 

22(it/24) 

.8550 

.6855 

1.926 

-.822 

-19.96 

1.926 

-.823 

-19.98 

.5000 

.4009 

2.434 

.368 

-24.49 

2.434 

.367 

-24.52 

24(  it/24) 

.8838 

.7086 

2.046 

-.871 

-24.62 

2.046 

-.871 

-24.63 

.5000 

.4009 

2.635 

.268 

-28.58 

2.635 

.266 

-28.59 

26(it/24) 

.9097 

.7293 

2. 157 

-.928 

-29.46 

2.158 

-.929 

-29.47 

.5000 

.4009 

2.829 

.153 

-32.89 

2.828 

.152 

-32.90 

28(it/24) 

.9326 

.7477 

2.261 

-.993 

-34.45 

2.261 

-.994 

-34.46 

.5000 

.4009 

3.013 

.024 

-37.40 

3.012 

.022 

-37.41 

30(it/24) 

.9524 

.7636 

2.356 

-1.064 

-39.56 

2.356 

-1.066 

-39.57 

.5000 

.4009 

3.186 

-.120 

-42.08 

3.185 

-.122 

-42.09 

32  ( it/24) 

.9690 

.7769 

2.443 

-1.143 

-44.76 

2.443 

-1.144 

-44.77 

.5000 

.4009 

3.346 

-.278 

-46.93 

3.346 

-.279 

-46.94 

34(it/24) 

.9822 

.7875 

2.521 

-1.228 

-50.02 

2.521 

-1.229 

-50.01 

.5000 

.4009 

3.492 

-.449 

-51.93 

3.492 

-.450 

-51.94 

36(it/24) 

.9919 

.7953 

2.590 

-1.318 

-55.27 

2.590 

-1.320 

-55.27 

.5000 

.4009 

3 1 623 

-.632 

-57.08 

3.623 

-.633 

-57.10 

38(it/24) 

.9979 

.8001 

2.650 

-1.414 

-60.49 

*  2.650 

-1.415 

-60.47 

.5000 

.4009 

3. ,736 

-.826 

-62.44 

3.736 

-.828 

-62.49 

40(it/24) 

1.0000 

.8018 

2.701 

-1.514 

-65.55 

2.701 

-1.516 

-65.51 

.5000 

.4009 

3.830 

-1.030 

-68.16 

3.829 

-1.034 

-68.37 

42(it/24) 

1.0000 

.8018 

2.743 

-1.619 

-70.32 

2.744 

-1.620 

-70.29 

.5079 

.4072 

3.901 

-1.242 

-74.80 

5.901 

-1.244 

-74.75 

44(it/24) 

1.0000 

.8018 

2.777 

-1.726 

-74.81 

2.777 

-1.727 

-74.78 

.5293 

.4243 

3.947 

-1.455 

-81.03 

3.945 

-1.456 

-80.91 

46(it/24) 

1.0000 

.8018 

2.803 

-1.836 

-78.97 

2.803 

-1.837 

-78.96 

.5599 

.4489 

3.969 

-1.661 

-86.33 

3.969 

-1.662 

-86.23 

48(it/24) 

1.0000 

.8018 

2.820 

-1.947 

-82.82 

2.821 

-1.948 

-82.80 

.5962 

.4780 

3.974 

-1.856 

-90.69 

3.975 

-1.856 

-90.59 

50(it/24) 

1.0000 

.8018 

2.831 

-2.059 

-86.28 

2.831 

-2.059 

-86.26 

.6354 

.5094 

3.966 

-2.038 

-94.16 

3.967 

-2.040 

-94.09 

52(it/24) 

1.0000 

.8018 

2.835 

-2.171 

-89.37 

2.836 

-2.172 

-89.34 

.6754 

.5415 

3.950 

-2.209 

-96.93 

3.951 

-2.210 

-96.88 

54(it/24) 

1.0000 

.8018 

2.834 

-2.283 

-92.07 

2.834 

-2.285 

-92.04 

.7149 

.5732 

3.927 

-2.369 

-99.11 

3.928 

-2.371 

-99.07 

S6(it/24) 

1.0000 

.8018 

2.827 

-2.396 

-94.40 

2.828 

-2.397 

-94.57 

.7531 

.6038 

3.900 

-2.520 

-100.83 

3.901 

-2.522 

-100.80 

58(it/24) 

1.0000 

.8018 

2.817 

-2.508 

-96.39 

2.817 

-2.509 

-96.36 

.7894 

.6329 

3.871 

-2.663 

-102.17 

3.872 

-2.665 

-102.14 

60(it/24) 

1.0000 

.8018 

2.803 

-2.619 

-98.05 

2.803 

-2.620 

-98.03 

.8235 

.6602 

3.841 

-2.799 

-103.19 

3.842 

-2.801 

-103.17 

62(it/24) 

1.0000 

.8018 

2.786 

-2.731 

-99.44 

2.786 

-2.732 

-99.42 

.8550 

.6855 

3.809 

-2.929 

-103.95 

3.810 

-2.931 

-103.94 

64(it/24) 

1.0000 

.8018 

2.766 

-2.841 

-100.56 

2.767 

-2.842 

-100.55 

.8838 

.7086 

3.777 

-3.055 

-104.49 

3.776 

-3.056 

-104.48 

66(it/24) 

1.0000 

.8018 

2.744 

-2.952 

-101.47 

2.745 

-2.953 

-101.45 

.9097 

.7293 

3.746 

-3.176 

-104.84 

3.746 

-3.176 

-104.83 

68(it/24) 

1.0000 

.8018 

2.721 

-3.062 

-102.18 

2.722 

-3.063 

-102.17 

.9326 

.7477 

3.714 

-3.294 

-105.03 

3.715 

-3.296 

-105.03 

70(n/24) 

1.0000 

.8018 

2.697 

-3.172 

-102.73 

2.698 

-3.173 

-102.72 

.9524 

.7636 

3.683 

-3.409 

-105.10 

3.684 

-3.411 

-105.09 

72(it/24) 

1.0000 

.8018 

2.672 

-3.281 

-103.14 

2.672 

-3.283 

-103.14 

.9690 

.7769 

3.653 

-3.522 

-105.05 

3.654 

-3.524 

-105.04 

74(it/24) 

1.0000 

.8018 

2.646 

-3.391 

-103.45 

2.647 

-3.392 

-103.44 

.9822 

.7875 

3.623 

-3.633 

-104.91 

3.624 

-3.635 

-104.91 

76(it/24) 

1.0000 

.8018 

2.620 

-3.500 

-103.66 

2.620 

-3.501 

-103.66 

.9919 

.7953 

3.594 

-3.744 

-104.71 

3.595 

-3.746 

-104.72 

78(it/24) 

1.0000 

.8018 

2.593 

-3.609 

-103.81 

2.594 

-3.611 

-103.79 

.9979 

.8001 

3.565 

-3.853 

-104.49 

3.566 

-3.855 

-104.49 

80(it/24) 

1.0000 

.8018 

2.566 

-3.719 

-103.90 

2.567 

-3.720 

-103.90 

1.0000 

.8018 

3.537 

-3.962 

-104.28 

3.538 

-3.964 

-104.30 

82  ( it/24) 

1.0000 

.8018 

2.539 

-3.828 

-103.97 

2.540 

-3.829 

-103.97 

1.0000 

.8018 

3.510 

-4.071 

-104.18 

3.510 

-4.073 

-104.19 

equation 


